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GEOMETRIC  DISCONTINUITIES  IN  ELASTOSTATICS1 

fay 

M.  K.  Kassir^  and  G.  C~  Sih^ 

Abe  tract.  Three-dimensional  elaetoetatic  problems  for  an  infinite 
solid  with  geometric  diecontlnuitiee  are  formulated  and  eolved  with  the  aid 
of  potential  functions.  In  the  problem  of  a  linearly  varying  preaeure 
specified  over  a  plane  region  bounded  by  an  ellipse,  uee  ia  made  of  the 
gravitational  potential  at  an  exterior  point  of  a  homogeneoue  elliptical 
diek.  The  problem  of  prescribing  displacements  on  the  eurfacee  of  disconti¬ 
nuity  ie  governed  by  the  Newtonian  potential  of  a  simple  layer  distributed 
over  a  diek  in  the  shape  of  the  region  of  discontinuity.  The  mass  density  of 
the  disk  is  proportional  to  the  prescribed  normal  displacement.  For  an 
"elliptically-ehaped"  region,  the  applicetion  of  the  symmetrical  form  of 
ellipsoidal  coordinates  leads  to  an  integral  equation  of  the  Abel  type  for 

the  potential  function.  It  ie  shown  that  if  the  displacements  normal  to  the 

2  2 
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elliptical  plane  are  given  by  (1  -  —g  -  Q  (x  ,y  ),  where  Qp(x  ,y  )  ie  a 

2  ga  b 

polynomial  of  degree  n  in  x  end  y  ,  then  the  corresponding  normal  etreec 

2  2 

acting  over  the  ellipse  ie  aleo  a  polynomial,  Pn(x  ,y  ),of  the  eeme  degree  in 
2  2 

x  and  y  .  In  the  caee  of  e  circular  region  of  discontinuity,  the  solution 
cen  be  carried  out  for  any  arbitrary  value  of  n.  The  reeulte  ere  useful  in  the 
prediction  of  the  stability  behavior  of  eleetic  eolide  containing  geometric 
discontinuities. 

*The  work  reported  in  this  paper  wae  carried  out  under  Contract  Nonr-610(06) 
with  the  Office  of  Naval  Research,  U.  5.  Navy. 
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1.  Introduction.  The  redistribution  of  stresses  due  to  the  presence  of 
geometric  discontinuities  in  elastic  aolida  has  been  the  subject  of  many 
past  discussions.  Among  the  three-dimensional  problems  considered  previously, 
the  majority  of  them  pertains  to  geometries  end  applied  loadings  that  possess 
exial  symmetry.  A  collection  of  previous  work  can  be  found  in  the  notes  of 
Sneddon  [lj  and  the  references  therein.  The  non-symmetric  problems  involving 
geometric  discontinuities  are  more  difficult  and  have  been  investigated  only 
in  a  limited  number  of  special  caaea. 

A  knowledge  of  the  clessicel  potential  theory  is  pertinent  to  the 
formulation  of  boundary  value  problems  with  planes  of  discontinuities  in  the 
three-dimensional  space.  For  definiteness  sake,  the  region  of  discontinuity 
will  be  assumed  to  take  the  form  of  a  plane  ellipse  described  by 


=  O 


(1.1) 


which  is  referred  to  a  set  of  cartesian  coordinates  x,  y,  and  z.  The  major 
and  minor  semi-axes  of  the  ellipse  are  denoted  by  a  and  b,  respectively.  If 
pressures,  p(x,y),  are  applied  symmetrically  to  the  upper  and  lower  sides  of 
the  elliptical  plane,  the  problem  reduces  to  the  determination  of  a  single 
harmonic  function  [2J  ,  f(x,y,z),  satisfying  the  boundary  conditions 


7ixt  0  /«  J  v>  A  *  * 


0 , 
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(~  4  1  >  I  *  =  0  ) 
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(1.2) 


(1.3) 


where  p.  is  the  sheer  modulus  of  the  elastic  solid.  Since  f(x,y,z)  is  harmonic, 


theae  conditions  mey  be  expreesed  in  the  equivalent  form 

,  (iti+p*1'}'0) 


(1.0 
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(1.5) 
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and  p*(x,y)  is  a  particular  solution  of  the  Poisson  squation  in  tun-dimensions 

2  jL  (1.6) 

The  problem  is  to  determine  the  potential  function  f(x,y,z)  such  that  it 

aatiefiee  (1.4)  and  (1.5). 

Green  and  Sneddon  £3]  have  examined  the  problem  of  a  flat  elliptical 
discontinuity  opened  by  a  uniform  pressure,  p(x,y)  =  constant.  In  part  of  the 
work  to  follou,  p(x,y)  ie  taken  to  be  a  linear  function  of  x  and/or  y,  i.e., 


-  C,  +  C2  %  +  C3  y  .  d.7) 

uhere  c  ,  (j*l,2,3)  are  constants  related  to  the  inteneity  of  the  applied 
preesure.  An  Bxact  eolution  ie  given  for  the  caee  of  Cj  =  c^  *  0  and  c^  = 
-Pq/2 jl .  The  eame  method  of  eolution  may  be  used  for  ci  *  c2  “  0  and  c^  * 
-pQ/2 jjL .  An  attempt  hae  aleo  been  made  to  find  f(x,y,z)  corresponding  to  a 
general  expreesion  of  p(x,y),  say 

m 

~  (  An  W  +  B*  ^  )  (1.8) 

ns  o 

The  end  reeulte  appear  to  be  rather  restricted  in  the  eenee  that  the  corffi- 
cienta  P|f  Bn  of  the  polynomial  p(x,y)  are  required  to  depend  upon  each 
other. 

An  alternative  formulation  of  the  problem  is  to  determine  the  distribu¬ 
tion  of  preesure  neceesery  to  preserve  certain  shape  of  the  eurfacee  of 
discontinuity.  Thie  can  be  accomplished  by  specifying  equal  and  opposite 
normal  displacements,  q(x,y),  on  the  region  given  by  (1.1).  For  thie  class 
of  problems,  the  boundary  conditions  are 
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where  0  ie  Poieeon's  ratio.  Equetione  (1.9)  and  (1.10)  show  that  f(x,y,z)  is 

equivalent  to  the  potential  of  a  eimple  layer  distributed  over  an  elliptical 

diek.  The  maee  deneity  of  the  diek  ie  proportional  to  q(x,y).  For  this  reason, 

eolutione  of  considerably  general  in  nature  may  be  found.  In  fact,  if  the 

o 

normal  displacements  are  given  by  a  polynomial  q(x.y)  of  degree  n  in  x  and 
2 

y  ,  then  the  integrals  representing  f(x,y,z)  cen  be  evaluated  in  cloeed  form. 
The  only  limitation  for  the  ellipticel  geometry  ie  thet  the  degree  of  q(x,y) 
must  be  equel  to  m+£y  m  being  a  positive  integer.  The  circular  rjgion  problem, 
however,  ie  free  from  euch  a  restriction. 

The  problem  of  specifying  pressures  and/or  displacements,  equal  in  magni¬ 
tude  but  opposite  in  sign,  to  the  plenee  of  a  flat  ellipee  may  also  be  formu¬ 
lated.  Kessir  and  Sih  (j»J  have  shown  that  the  the  ekew-symmetric  problem 
requires  the  knowledge  of  two  hermonic  functions  g(x,y,z)  and  h(x,y,z).  The 
simple  ceee  of  en  elllpticelly-shaped  plane  of  discontinuity,  whose  eurfecee 
ere  subjected  to  uniform  sheer,  wee  diecuseed  in  [4j  •  Without  further 
comments,  it  ie  obvioue  that  a  generel  treetment  of  the  skew-symmetric  problem 
followe  immediately  from  the  work  presented  in  thie  paper  end  in  [ 4 J  . 

2.  Equetione  in  Elestoetatlca.  Consider  the  problem  of  finding  stresses 
end  di8plecemente  in  an  infinite  elastic  eolid  whoee  continuity  is  interrupted 
by  a  void  in  the  ehape  of  a  plene  ellipee.  The  epplied  pressures  or  displace¬ 
ments  on  the  upper  eurfece  of  thie  ellipee  are  equal  to  those  on  the  lower 
eurfece.  Hence,  the  problem  ie  said  to  be  symmetric  with  respect  to  the 
ellipticel  region  located,  eey  at  z  ■  o  or  in  the  xy-plane.  Outeide  of  this 

region,  the  eheer  etreeeee  T  end  T  must  vanish  et  z  »  o.  These  conditions 

x  z  yz 

mey  be  satisfied  by  expressing  the  certeeien  components  u,  v,  and  w  of  the 
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displacement  vector  in  terms  of  a  single  harmonic  function  f(X),y,z)  as 


+  , 

(2.1) 

7i + }  ' 

(2.2) 

)  -f  -y.  Ol.  . 

(2.3) 

here  since  they  have 

already  been  given  elsewhere  [2] 

V 


In  the  usual  way,  the  stress  components  l~xzt  7^,  and  ^zz  arB  9lven  by 


7^  -  2/lfr  - „  y 

r» 
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Tul  „  2 
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(2.4) 

(2.5) 

(2.6) 


W 

-  ^£n 

For  the  mere  purpose  of  determining  the  potential  function  f(x,y,z)5  thB 

remaining  stress  components  ^  ,  6"  .  Bnd  T  are  not  needed. 

xx  yy  xy 

3.  Linearly  Varying  Pressure,  lei  the  two  sides  of  the  elliptical  plane 
of  discontinuity  be  subjected  to  a  pressure  thBt  vBrie3  linearly  with  the 
coordinate  x  over  the  region  in  (1.1).  The  disturbance  of*  the  discontinuity 
diminishes  at  large  distance  awBy  from  the  origin.  On  the  plane  z  =  o,  (2.3) 
and  (2.6)  yield 

is*-  v  a  bl  '  /  J 


*3 

if  -  °. 


2/1 


(3.1) 

(3.2) 


It  is  assumed  here  that  the  second  and  third  derivatives  of  f(x,y,z)  with 
respect  to  z  are  bounded  in  the  limit  as  z  approaches  zero.  These  require- 
ments  will  be  met  in  the  subsequent  analysis. 
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The  problem  can  be  resdily  solved  upon  defining  a  new  harmonic  function 
<f>(xtytz)  related  to  f(xyyyz)  by 


where 


H 
n 

v fa-Vi)  =o  ,  -  o 


#  *A  -  yl£  , 

*  9?  *7# 


(3.3) 


and  V  ie  the  Laplacian  operator  in  three  dimensions.  From  (3.3) y  it  can  tie 
shown  that 


w  =  +  £  >  % 


f  }  ryi 


£f  _  3L, 

9f  " 

The  limite  of  the  integrals  in  (3.4)  snd  (3.5)  sre  introduced  in  euch  a  way 


(3.4) 

(3.5) 


that  the  stresses  snd  displacements  vanish  ss  z  goes  to  infinity.  In  terms  of 
the  function  ^(x,yyz),  the  displacements  became 

oo  oo 

VL  =  0-20)  +{}%  dj]  +  6) 

r  -  W)  $  0  f  ^  O  (3>7) 

V  -  - j  -;g]  • 

Using  (3.3)-(3.5) ,  the  stress  components  sre 


T 


1} 


(3.9) 

(3.10) 

(3.11) 


o 


Thus,  ths  boundary  conditions  (3.1)  and  (3.2)  taka  the  foro 


* 


1  9% 


(fi  +  <  1  >  }"  o) 


(3.12) 


(3.13) 


i  ov. -  * 

Upon  introducing  a  system  of  ellipsoidal  coordinates,  the  position  of  a 
point  can  bs  determined  by  the  three  parameters  E, 7  ,  C,  ths  limits  of  varia¬ 
tion  of  which  are 


<*>>■£  io  ^7  >  -bl  >  Z  ^  -a*\  (3.i4) 

The  cartesian  coordinates  x,y,z  may  be  expressed  in  terme  of  theee  parameters 
by  ths  relations 


a  (a-  bl)  %  =  (a  +  i)(a  +  rj)(a+C) 

(3.15) 

i  (b-tC)  y  s  (b  *  £  )(fc  +  (  b  +  O  / 

(3.16) 

«V;‘  -  £75. 

(3.17) 

2  2  2  2 

In  the  plane  z  -  o,  the  inside  of  the  ellipse  x  /a  +y  /b  =  1  may  be 
distinguished  from  ths  outside  by  setting  E  =  °  and  *  o  ,  respectively. 

The  boundary  of  the  ellipee  is  identified  by  ^  ^  =  o, 

Since  the  function  ^(x, y,z)  is  equivalent  to  the  gravitetional  potential 
of  e  homogeneous  elliptical  disk  at  an  exterior  point  x,y,z,  the  solution  ie 
given  by  [5] 

<x> 


where 


(3.18) 


Q( 5)  «  S(«V S)(b\s)t 


(3.19) 


Inserting  (3.18)  into  (3.12)  and  (3.13),  ths  constant  A  is  found  to  be 


h 


4!1  KMU-  Ed) 


(3.20) 


In  (3.20),  K(k)  and  E(k)  are  Legendre's  complete  elliptic  integrels  of  the 
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first  and  second  kind,  respectively,  for  the  modulus 

*l-  /-  (jfj\ 

2  2  jl 

and  k'  =  1  -  k  .  Once  9<x,y,z)  is  known,  the  stresses  and  displacements  at 
any  point  of  the  aolid  can  be  calculated  in  a  straightforward  manner. 

For  instance,  the  displacements  u,  v,  and  u  in  the  plane  z  =  o  may  be 
obtained  readily  from  (3.6)-(3.8)  which  reduce  to 

ti  -  (1-10)  fg(**+gf), 


v 


-  20-°)  'X,  §£ 
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(3.21) 

(3.22) 

(3.23) 


in  which 


OO 


(3.24) 


(3.25) 


tcxyt)  «  jT  >  , 

Introducing  the  contraction 

.Vsl  i  _  _  _2l  _  i!  , 

feVs  S 

and  substituting  (3.18)  into  (3.24),  the  function  ^(.x,ytz)  may  be  represented 
by  the  integral 

Vh-*»--}L>Ua»dkh 


_  A 


(3.26) 


It  follows  that 


4  4 [m>]  1  Mk 


»y 

i 


06 


**  ’  *  4  ~ k  4  f^N/OS) 


Hl 


06 


-/)»*  /  r  5 4s 

9  frTi  e  \/L*i 


(3.27) 


(3.28) 


3#^  '  "T  (0l+s)(bz+s)jQCs) 

The  results  of  putting  (3.18),  (3.27),  and  (3.28)  into  (3.21)-(3.23)  are 


(3.30) 

(3.31) 


oo  .<x>  ^  06  d  S 

U  '  (,~l0)  lt  [3*i6  +  ficalsxfetffc 

oC 

V-d-ZibAa.1**  \  ,  lte.  tti  , 

'  J5  (it%sXfcVs;>/Qs; 

if-  -Z(l-\»/)x  Um.  [%  f  ^=]  * 

The  integrals  in  (3.29)-(3.31)  may  be  evaluated  in  closed  form.  They  are 

04 

06 

f  of£  ^  [  r,  ,  u>z  Mi*  CJIU'  1 

4  (b*+sjJQ&)  "  a’-*.1-*'1  ;  u~  *  4n  il  J  ' 

and  others  of  the  same  type.  The  quantity  E(u)  is 

U 

E(u)  **  J  dn2t  dt  / 


29) 


(3.32) 

(3.33) 


(3.3/*) 


and  snu,  cnu,  dnu  are  the  Jacobian  elliptic  functions  whose  variable  u  should 


be  distinguished  from  the  x-component  of  the  displacement  vector.  Here  u  is 
related  to  the  ellipsoidal  coordinate  ^  by 


(X  an 


V/  snu 


a2(Snu 


-0. 


(3.35) 


After  some  manipulations,  the  displacements  for  z  *  o  are 

U  ~  (-$rr{xl[(,4-%r)u  -z(l+jp)E(u)+  Snu  cnu,  dn  itj 
-  yl[j?U-^rfcrE(u)  +  ^snitcmflnu.]  +  al[E(u)-u]J  >  '.3.36) 


^  t6  “1 

'j£z  cJu  IL  J 


/ 


(3.37) 
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(3.3B) 


On  the  plane  surface  of  the  discontinuity,  i.e.,  for  £  *  o,  (3.36)-(3.38) 
may  be  simplified  since 

fiM-Eflb  ,u~K(i)  ,  — o- 

^  1  2  2  2  2 

Hence,  the  deformed  shape  of  the  ellipse  x^/a^+y  /b^  =  1  is  given  by 


(H«o"  { ^[(z+41)  k(t)  -2(1+4}) E(i)] 

-yl[*K(t)-U+ -&)£«)]  +a*'kt[E(t)-K(k))} ' 

[('+&)bu -z  m] , 


(3.39) 


(3.40) 


_ ifMM  i » 


(3.41) 


5’«  ab  b' 

The  constant  A  hss  alresdy  been  determined  in  (3.20). 

Of  interest  is  the  intensification  of  the  normal  stress  6*22  in  the 
neighborhood  of  the  geometric  discontinuity.  It  is  the  local  1  ovation  of 
6*  in  the  vicinity  of  the  edge  of  the  discontinuity  (  %  -►t*,  0  -*a)  thBt 
controls  the  stability  behavior  CO  of  the  plane  elliptical  cavity.  From 
(3.11) 

%  “  ~xiy)  >  ».«> 

Inserting  (3.18)  into  (3.42)  and  integrating,  the  obtained  result  gives 


(3.43) 
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The  expression  (3.43)  is  valid  only  in  the  region  outside  of  the  ellipse 
2  2  2^ 

x  /a %y  /b4"  *  1,  To  find  0”  near  the  periphery  of  the  elliptical  boundary, 


& 


the  asymptotic  value  of  2*  is  required.  By  taking  a  vector  f  normal  to  the 
curve  that  defines  the  ellipse  with  parametric  equations 

%  »  a.  sin(j>  ,  y  *  b  C03<f> 


(3.<44) 


the  expansion 


£a.b  p _ 

5  “  (a*sin*<f>  +  blcc6%)'M  +  ""  (3.45) 

is  obtained  for  which  terms  of  order  higher  than  p  have  been  neglected 

since  p  ,  the  magnitude  of  p  ,  has  been  assumed  to  be  small  in  comparison 

uith  either  a  or  b.  In  addition,  it  can  be  shown  that 


(3.46) 


(3.47) 


Z  =  — (<LStri'4>  4-  bl CoSl4>) 

The  limiting  forms  of  £  and  £  lead  to 

(6>4-°  4  cos^ 4  •  •  •  > 

Denoting  the  coefficient  of  l/\f 2p  by  k1  and  using  (3.20)  yield 

This  is  known  as  the  three-dimensional  "stress-intensity  factor"  [*»3  the 
critical  value  of  which  governs  the  onset  of  unstable  motion  of  the  geometric 
discontinuity  under  consideration. 

2  2  2  2 

The  applied  pressure  on  the  ellipse  x  /a  +y  /b  =1  may  also  depend  on 
^2n  y2n^  unere  n  is  a  positive  integer.  For  n  =  1,  the  function 


oo 


-  xl 


ds 


*5  aj?&) 

will  satisfy  the  boundary  conditions 

^  E(k)  -  [(jt-O £<*)+«£/]  % 


(3.49) 


(3.50) 


■=■  11  «■ 


(3.51) 


-i- 


-  jj-n  [(  ^&c)EG)  "K(*)J  yl|  >  £  -0 


3}  -o  ,  7-0 

However,  the  result  baaed  on  (3.49)  will  obviously  be  United  to  the  type  of 
pressure  distribution  as  indicated  in  (3.50).  Similar  observations  may  be 
made  for  n  greater  than  one. 


4.  Specification  of  Displacements.  As  mentioned  earlier,  when  the  region 
of  discontinuity  undergoes  certain  normal  displacements,  q(x,y),  that  are 
preacribed,  the  problem  can  in  general  be  solved  for  the  diatribution  of  the 
normal  pressure  which  is  initially  unknown.  The  conditions  (1.9)  end  (1.10) 
suggest  taking  the  function  f(x,y,z)  as  the  Newtonian  potential  of  a  simple 
layer  of  intensity  q(x,y)  distributed  over  the  plane  region  Q  in  the  shape 
of  the  ellipse  x2/a2+y2/b2  *  1.  The  potential  is 


•#*»>  -  «b,  £f  ^  • 

where  R  denotes 

R  -  [(v-Xf+ty-Y)'-*  f]  . 

The  function  f(x,y,z)  ia  harmonic  and  ia  continuous  in  the  entire  three- 
dimensional  region  excluding  fl  .  Furthermore,  it  vanishes  for  distances 
sufficiently  far  away  from  since 


4--k[w- s^lVwxdy]  ,  « 

The  normal  derivative  of  f(x,y,z)  ia  diacontinuous  for  the  transition  from 
one  side  of  .0.  ,  say  z  =  o+,  to  the  other,  z  &  o”,  i.e., 


(4.2) 

(4.3) 
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Hence,  (1.9)  and  (1.10)  are  satisfied  automatically  when  f(x,y,z)  is  cons¬ 
tructed  from  the  potential  of  a  simple  layer. 

For  ths  present  problem,  it  is  more  convenient  to  employ  the  symmetrical 
form  of  ellipsoidal  coordinates  5,  tj  ,  £  as  given  in  (3.10-0.17)  • 

Proposed  as  a  possible  solution  is  the  function  [£j 

oo 

W  *  6jfS_ 

(4.4) 


-  1 


Ths  variable  CO  is  defined  by  (3.25)  or  the  equivslsn*  ^orm 


(S-Z)(S-0)(S-Z) 

Q(s) 

and  A  is  a  twice  diffsrentiable  function  in  the  interval  (o,l)  with  finite 


Ms) 


one-sided  derivatives  st  the  boundary  points  of  the  interval.  Differentiating 
(4.4)  with  respect  to  x  renders 


Hi  _  fm?)  ds  4  Mo)  , 
ti.  f  is_  .  u.  £*  2i  +  >&) 

?**■  4  «x  JQfi;  3%  u 


ii&) 


(4.5) 


(4.6) 


in  which 


oo  . 

r 


satisfies  the  harmonic  equation  *  0.  The  remaining  expressions  of 

2  2  2  2 

*df/4dy,  d  f/'Oy  ,  and  Vf/dz,  7l  f/*dz  ,  may  be  obtained  simply  by  permu¬ 
tation  of  tha  variable  x  in  (4.5),  (4.6)  to  y  and  z,  respectively.  As  s 
result,  the  Laplacian  of  f(x,y,z)  can  be  writted  as 

oo 


vV*' w  -  {  4  ^  m> 

with  the  knowledge  that  [7] 

^  ?I  +  a  ;I  ±  UL7£ 

Now,  computing  the  derivatives 


(4.7) 


I  . 


13  - 


9% 


z% 


al+s  ^'(a))  * 


(4.8) 


ax*  fl*+s  * 

and  etc.,  the  integrand,  V^A(CO),  in  (4.7)  becomes 

By  means  of  (4.7)  and  (4.8),  f(x,y,z)  is  indeed  a  solution  of  the  equation 

There  remains  the  determination  of  A  (a))  which  may  be  evaluated  from 
the  boundary  conditions  (4.2)  and  (4.4).  First  of  all,  replacing  x  by  z  in 
(4.5)  and  using  the  relation 

Ik  _  itQgl 

-L 


the  result  is 


(4.9) 


Uhile  (4.9)  satisfies  (4.3)  for  ^=0,  it  reduces  to  an  integral  equation  of 
the  Abel  type  for  2;  *  o  as  dictated  by  (4.2).  This  will  be  shown  in  the  next 
section. 

5.  The  Abel  Integral.  Applying  the  boundary  condition  (4.2)  for  points 

2  2  2  2  4> 

of  xp  y  inside  x  /a  +y  /b  =  1  and  z  -»■  0  ,  (4.9)  may  be  written  as 

X 

°kX(\y  -  tim,  [($qefj7&» 

tLfL- ill  .  A  4  »  r 


c/s  1  a.b  7i(o) 

'sjQfe)  J 


(5.1) 


Wfa)  z+o  u  •  £ 

The  reduction  of  (5.1)  to  Abel’s  equation  may  be  accomplished  by  changing  the 
variable  of  integration  from  s  to  o(  *  1-5/s  and  by  keeping  in  mind  that  60 
is  a  function  of  s  with  the  limiting  form 

7$  ■  os 


axs) 


where 


*7  C  -►  c?  bx  Z  ,  and.  Z  *  /“  ^  * 


(5.2) 


(5.3) 


Under  these  considerations,  (5.1)  reduces  to 

4m2  ' z*/.  4  ^ 

An  additional  change  of  variable  to  ^  =  o(  Z  gives  the  Abel  integral  equation 

[6]  ,  2  /#  , 

£jj(Z)  f  d?  . 

J0  Tz^r  C5.<.) 

In  order  to  solve  for  ^.(0)),  both  aides  of  (5.4)  are  integrated  from  o  to  60 
by  dZ/yTu)  -  Z  ’  •  Hence,  (5.4)  becomes 


(W 


1 


Xo)  dZ 

Vi 


f  jl.  r  ‘t&jz ' 

IJZFZ  J„  tWw5  "'o 


a) 


(5.5) 


K  LUt-s))J0  sTZTz 

The  integrala  on  the  left  hand  side  of  (5.5)  may  be  evaluated  uithout  diffi¬ 
culty  the  first  of  which  ia 

U)  X 

f  «  f  Tfexto  =  TrXp) 

Jo  Z  Jw-Z  7o 

The  order  of  integration  of  the  second  integral 


(5.6) 


f  _iz_  f  ’aWe 

-L  \Jfj-r7.  Jo  vTZ -  ^ 


(5.7) 


o  CJ-*7,  ^  o 

may  be  interchanged  since  the  integrands  possess  weak  singularities.  Applica¬ 
tion  of  the  Dirichlet  formula  [7]  gives 

6)  0)  uj 

dZ 


i>^i 


nrf  ?i'c$)Jp  «  t[?)(co)-7)(o)] 


(5.B) 


Substituting  both  (5.6)  and  (5.B)  into  (5.5)  yields  the  solution  of  (5.4)  as 

^  =  i^hjl  •  (5-9) 

Thus,  (5.9)  solves  the  problem  of  prescribing  normal  displacements  q(Z)  on 
an  elliptlcally-shaped  plane  of  discontinuity. 

6.  Theorem.  On  the  basis  of  the  results  obtained  in  *he  previous  section, 
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* 


the  following  theorem  may  be  established: 


"Let  the  normal  displacements,  q(x,y),  on  the  plane  surfaces  of  the 

2  2 

ellipse  x2/a2+y2/b2  =  1  be  given  by  (1  -  QnCx28y2)  v  where  Q^x^y2) 

2  2  a  ^ 

is  a  polynomial  of  degree  n  in  x  ,  y  .  Then,  the  normal  pressure  acting  over 

2  2  2  2 

the  ellipse  is  also  a  polynomial,  P(.(x  ,y  ),  of  the  same  degree  in  x  ,y 


Now,  suppose  that 


®  Z  Z  Z  . 


(6.1) 


2  2, 


It  suffices  to  prove  the  theorem  by  taking  Q  (x  ,y  )  in  the  form  of  (6.1). 

The  coefficients  Cj  ( j=o,l, • • • ,n)  are  constants  and  Z  is  given  by  (5.2). 
Putting  (6.1)  into  (5.9)  and  carrying  out  the  integration  yields 

?(co)  =  mjf  f:0  -T;if f  Ci  ^  *  <*■» 

in  which  P (n)  is  ths  customary  Gamma  function.  Before  the  normsl  pressure 
6 z2  csn  be  calculated,  the  harmonic  function  f(x,y,z)  must  be  obtained  from 


(4.4): 


ab 


Cj  '  C6.3) 


In  the  plane  z  *  a,  (2.6)  simplifies  to 


(6.4) 


Adopting  the  notstion 


7  = 


0  ,  <l>  }*°) 
6.  (-5-?>»^-) 


(6.5) 


and  using  (6.3),  it  is  found  thst 

ft 


i  j(^4 (6<6) 


where 


i^)*N},0-  '-  -  -fc ' 


The  integrands  containing  (ft)  (a)]n,  (n**l, 2,  •••)>may  be  expanded  in  terms  of 


2  2 

x  ,y  to  give 


»  f\,  *  A,  -  -  -  -  P.ftC  ,  (6.7) 

in  which  Aj  C j*o, 1,2, * • *fn)  are 

Ro  -  6.J  («4s  +  055 )  jfit)  ’ 

“  6,X  (a+s)(b%s)lQ<s)  *  &tIr  (£+* ^sjGps)  ’ 

#*•  *  6'/  GVs HbU)Wt)  *  fitX  ’ 

and  so  on.  The  constants  0q,  •••,  09  are  related  to  Cj  (j=o,l,2)  as  follows: 


e>.-^/L j \Hhkic.  s, rM±h-c-  6,-36, 

0  O-otfr  poji  W  '  '  MiWfz,  (/->)[  i  >  2  ' 

_  2  2 
From  (6.7),  (  02Z)2_Q  is  seen  to  be  a  polynomial  of  degree  n  in  x  ,y  whose 


coefficients  depend  on 

do  of 


©o 


f 


dS  f  dS _  I  _ , 

(qVs)WH/Q(sj  ^  £+$)M\]Q  &)  JQfe) 


dS 


ds 


0,  (;-•  - 


These  integrals  can  be  reduced  to  the  complete  elliptic  integrals  of  the 
first  and  second  kind  with  modulus  k  =  1  -  (b/a)  for  £  *  o. 

The  aforementioned  theorem  also  applies  to  the  skew-symmetric  problem 
of  specifying  displacements  u ,  v  an  the  elliptical  plane  of  discontinuity. 
Thia  ia  mainly  because  the  boundary  conditions  on  u,  v  are  the  same  as  those 
described  by  (4.2)  and  (4.3)  since  Kassir  and  Sih  [4]  have  already  shown  that 
fur  the  skew-symmetric  caae  the  displacements  are 
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inhere 


U  =  -  ZCl-O)  ||-  •+  f  ^  > 

V--ZW 

vr~  -0-24)6.  +  »>!, 

*  ?y 

r  tl  +  &  , 

w-  -  fx 


v'4(z.^)  -  o. 


For  this  problem,  the  snear  stresses  T 9  T  _  will  be  polynomials  of  degree 

xz  yz 

2  2 

n  in  x  ,  y  if  the  displacements  u,  v  are  represented  by  polynomiale  of  the 
form  (6.1). 

In  the  two-dimensional  case,  Sneddon  [l]  hae  observed  that  the  preecribed 
ehepe  of  e  slit  in  an  isotropic  medium  end  the  resulting  eurface  pressure  may 
be  represented  by  polynomials  of  the  seme  degree  in  one  varieble.  Galin  C&3 
has  considered  tr.e  three-dimensional  problem  of  en  ellipticel  punch  preesing 
against  a  semi-infinite  elastic  solid  whose  surface  outBide  thB  '>rea  of 


contect  is  free  from  trections.  The  frictional  forces  ecting  between  the  punch 
end  the  semi-infinite  body  are  neglected.  He  showed  that  if  the  base  profile 
of  the  punch  is  given  by  e  polynomial  of  degree  n  in  x,  y,  then  the  pressure 
ecting  over  the  punch  is  enother  polynomial.  Gelin'e  proof  ie  based  upon  the 

5 

properties  of  Lame  triple  product,  which  ie  different  from  that  given  in  the 
present  paper. 


7,  The  Symmetric  Problem.  A  sufficiently  generel  description  of  the 
shape  of  the  elliptical  plene  of  discontinuity  is 

C7.1) 

where  qQ  denotes  the  amplitude  of  the  normal  displacement.  The  restriction 
on  the  exponent  n  becomes  apparent  when  the  functions 
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—■*»**-  ***%. 


(7.2) 


and 


oo 

M%m) _ iW*2L.[  [o&f^jdL  ,  o-i  , 

0  iuh>) Jt  ri*+£y-  *■  J  >/35)  ' 


(7.3) 


are  derived  from  (5.9)  and  (4.0,  respectively.  The  normal  pressure  distri¬ 
bution  required  to  support  the  shape  (7.1)  can  be  obtained  from  (6.4)  and 
(7.3).  For  the  range  of  n  between  1/2  and  3/2,  6"  in  the  plane  2  *  o  is 

4t 


ffWw  -  » si  *  t  &U  ‘  ."•« 


V*  ‘'Z  V  •»  ■  v  ■  2 

where  I,,  I„  represent  the  integrals 

■"rW  .  x 


-  _  f  f&w 
■' '  i  r«i+ 


S) 


J 


C&W  1  </s 

(fc‘+s;  jOte) 


w 

and  T  is  defined  by  (6*5).  If  n  ^3/2,  the  expression  for  (Qzz\-0  ia 


less  complicated: 


AUU w  wumpA^uu  Ubu  a  ^ 

^-r  -  (I  <ss 4  M  fijl 

-  (u-oj  [(&)  +  (j£s)1]  I  • 


(7.5) 


In  general,  the  integrals  appearing  in  (7.4)  and  (7.5)  can  be  evaluated  for 
n  *  m  +  %,  where  m  -  o,l„2,  .In  the  degenerate  case  of  a  circular  region 

of  discontinuity,  solutions  may  bs  found  for  any  value  of  n. 


6 ,  Elliptical  Region  for  n  a  1/2  and  n  =  3/2.  A  special  case  of  (7.4)  is 

* 

n  a  1/2,  which  corresponds  to  the  problem  of  a  plane  elliptical  cut  opened 
by  a  normal  pressure,  ( 02J)  =  constant.  The  opening  of  the  cut  is  ellip¬ 
soidal  given  bv  ± 

This  problem  was  formulated  by  Green  ans  Sneddon  [3  ]  .  In  a  later  paper, 
Kassir  and  Sih  [4]  pointea  out  that  the  stresses  at  the  border  of  ths  cut 
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are  singular  of  the  order  of  1/f^* •  The  distance  ^  is  measured  in  a  plane 

2  2  2  2 

perpendicular  to  the  boundary  of  the  ellipse  x  /a  +y  /b  =1. 

Therefore,  it  is  natural  to  inquire  whether  cuts  of  a  shape  other  than 
ellipsoidal  are  possible  and,  if  they  exist,  to  determine  the  pressure  neceseary 
to  preserve  their  shape.  Consider  the  case 

3/z. 


=  f.  Z  . 

For  n  =  3/2,  (7.5)  specializes  to  ^ 

~  +  -™llr 


(8.1) 


~  z  c/s 

^  1  (feVsjVQs;  ^  ^ '  i  (<hs)(lte){Q($) 


]■ 


(8.2) 


The  integrals  in  (8.2)  may  be  evaluated  in  a  manner  similar  to  those  shown  in 
(3.32)  and  (3.33). 

The  normal  pressure  inside  of  the  elliptical  region,  T *  o  ,  is 


(8.3) 


When  a  *  b  and  E(k)  *  H(k)  =  'T/ 2 ,  (8.3)  gives  the  internal  pressure  for  a 
penny-shaped  cut  of  radius  s,  i.e., 


(8.4) 


The  sign  of  (  ^22)z=a  changes  from  negative  to  positive  st  the  ratio  of 
(r/a)  £*  CL6.  This  means  that  both  compressive  snd  tensile  stresses  must  be 
applied  to  the  surfaces  of  the  cut  in  order  to  produce  the  shape  (8.1). 

Outside  of  the  region  of  the  ellipse,  (  ^2Z^Z_0  rosy  be  found  by  setting 
*  £  in  (8.2) : 


(%V.' 


Snu  Cnu, 

Jnu. 


(8.5) 
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.  .  .»»,  m_.  .. _  .  ^  ..  .»tSi3S&fc 

jm 

’&&«+ o-*#)Bto  +  sag&^^-oj 

In  contrast  to  the  singular  solution  for  n  &  1/2,  the  stresses  remain  finite 
on  the  boundary  of  the  elliptical  cut.  In  fact,  in  the  limit  aa  %  -•*  o,  (8.5) 

JL 

is  of  the  order  of  ^  .  The  same  behavior  is  observed  for  the  particular  case 
of  a  *  b  since  (8.5)  reduces  to 


3, 


r>a 


(8.6) 


where  (  (T  )z  goes  to  infinity  ae  0(l/r  ).  When  r  approaches  a,  the  circular 
boundary,  the  normal  pressure 


(6,,)  -  1&2L  (Is)  r-a.  (8  7) 

#7-0  2o-v) v  a '  »  18  7) 

is  found  to  be  non-singular. 

It  should  be  mentioned  that  (9.4)  and  (8.6)  may  be  obtained  directly 
from  the  method  of  Hankel  transform  for  solving  axially  symmetric  problems  [9J  . 


9.  Circular  Region  for  Arbitrary  n.  The  theory  of  potential  functions  may 
also  be  applied  expediently  to  a  class  of  problems  involving  penny-shaped 
planes  of  discontinuities.  For  a  =  b,  (7.4)  can  be  solved  in  general  without 
imposing  restrictions  on  the  exponent  n  as  1^  =  I2  =  I  and 


The  integral  I  is 

1 


OO 

f  [(S-DCS-/))]  '  *  a's 
\  S*  0Vs/+T 


and  the  symbol  stands  for 


r 


rco. 
a  . 


(9.1) 


(9.2) 


(9.3) 


M 


i 


T 

1 


a 

i 


w 

A 

I 

w! 

J 

ai 


Ii 

ti 


-  21  - 


Making  the  substitutions 

r=  ia  > 


e~= 


s -1 

in  (9.2) ,  the  integral  I  becomes 

r  ?j2"’  f 

J  /-  _»  .a  ..  * 


(9.4) 


CaV*r*  Jo  Ctr+ifa-fisf** 

The  determination  of  (  ^ZZ^Z_Q  for  r<®  ®nd  r>  a  will  be  carried  out  separately, 

2  2  2 

Inside  of  the  region  x  +y  =  a  and  fc,r  z  »  of  the  ellipsoidal  coordinates 
%  ,  take  the  values  x  , 

5-o,  (\  =  r  -  a  , 

2  2  2 

and  €  =  k  3=  (r/a)  .  Hence,  (9.4)  reduces  to  the  standard  Legendre  elliptic 


integral 


where 


CD. 


*JlCr 

a3  J°  (i-fo1)"**  * 


I 


(9.5) 


4T-  1-  C£)‘. 

A  further  substitution  of  =  snu  leads  to 

wi..- 


(9.6) 


In  (9.6),  represents  the  integral 


i  f?L 

Jtn  *  J  W  a)  eta  , 

whose  values  for  m  =  o,l,2,  •••,  are  given  by  [lo] 

Jo  “  J  J"|  =  <  Jt" 


IT  T  _  £(&L  T 

/z  ^  a 


? 


(9.7) 


(9.8) 


and  the  recurrence  formulas 


J. 


I 


X 


zm+z  (2 

I 


[zm(£-'ti)Jvn+  0-2m)J^z  ]  , 


[(2fll+i)(Z-  it)  lint,  ~  Jam- 1  3 


(9.9) 


(9.1Q) 


zrn+z  zCmH)^7- 

In  order  to  complete  the  calculation  of  (  ^zz^2_0  for  £  =  o ,  it  is  necessary 


to  find  from  (9.2)  the  quantity 


rr«fc*»5uL-*r* 


(9.11) 


Some  lengthy  algebra  gives 


(?1)  JT^JS-  h„n)f  ^ 1 -$*■<&  1 

Vs-.  a*  L  **  l  jr±Zr  i  ™»>J  ~7T^Fl  J 


-_j£/ *2. 

a5-*1  (  ■*“ 


(HEV)'1  %  «-  *V/ 

>(^)[  J.-gh6  4* + *'24.j} 


(9.12) 


uihere  J  is  given  bv  (9.7).  Using  (9.1).  the  internal  pressure,  (  (T__)_  „  , 
n)  zz  zsd 


in  the  region  r<a  takes  the  *orm 


(9.13) 


(6ij  =  _  yj (n*2) Ban ci») f r_  (14^ ) j 

33  **°  u^JWW^fy«- ''**'*•  (  <#*  Jz 

+  (34-U-n.  -f  it";  J2fl«  -  (2n+i)  ^  ]  , 

which  is  valid  for  arbitrary  value  of  n.  Three  special  cases  of  (9.13)  will 
be  considered  by  application  of  (9.8)=>(9.10) . 

Casa  9.1(a) >•  The  necessary  normal  pressure  to  maintain  the  shape 


/>  /  ,  /  r  »2  1  t 

0 

may  be  computeo  from  (9*13)  by  setting  r.  =  1/2  as 


‘  +  (5+^  J*  -  ]  C9.U) 

>  r<4 

This  i8  in  agreement  with  Sneddon’s  answer  [9  J  • 

Case  9.2(a) .  If  n  =  1,  (9.13)  simplifies  to 

«H..-  -  zSftMtt-  i  *  f*  «  j.  -  mv,  ] 

At  first  sight,  (9,15)  appears  to  be  different  from  the  solution  published 


(9.15) 


so*. 
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by  Sneddon  [9]  for  the  same  problem,  i.e., 

<•*»&]  ■  <»•■» 


in  iiihich 


4C' 


Q-ar 
(Q+r) 1 


However,  by  means  of  the  Landen  transformation,  the  complete  elliptic  integrals 
in  (9,16)  can  De  written  as 


(1  +?)£(**>-  *  E(g) -[/-&]  k®  , 
KM  -  K(£f. 


(9.17) 


(9-IS> 

The  identities  (9.17),  (9.18)  show  that  (9.16)  and  (9.15)  are  indeed  equivalent, 
Case  9.3(a) .  In  the  case  of  n  *  3/2,  (8.4)  is  recovered  from  (9.13) 

since 


(9.18) 


£ 


In  the  same  fashion,  (6*  )  outside  of  the  circular  plane  of  discon- 

72  2=0 

tinuity  may  be  obtained  by  letting 


and 


«v.  - 


Z  'k 


tin 


( 


©*77- e1  Jr 
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which  is  found  from  (9.4)  for  ^  15  0  ,  €  =■  k  =  (a/r),  and  r  -  a  .  The 

substitution  (T  *  snu  allows  (9.20)  to  be 


(9.19) 


(9.20) 
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(II  =  £i_ 
h°  r**,1 


hUf  » 


(9.21) 


where 


m 


J 


(dnuf 


n+z 


(9.22) 


(9.23) 


The  integrals  (9.22)  can  he  expressed  in  terms  of  the  complete  elliptic 
integrals  and  elementary  functions  in  the  following  manner  [loj  : 

kA  dfl'(4)-sitl(m^)  k  a  E(i)-^lKlit)  hj  M, 

4-k'  ’  *“  4l4l  '  ' 3 

The  recurrence  relationships  are 
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In  addition,  the  integral  (9.23)  admits  the  relations 
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If  the  exponent  n  -  m/2,  m  Deing  an  integer,  tnen  (9.20)  is  reducible  to 
elementary  ^unctions.  Finally,  the  normal  pressure  for  r^a  is  deduced: 

(SjJ  c  -  ty(n+2),'top±/}9\j'2n  f z  ,i  _  J,  ♦  „ViXJ  ) 
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The  valuee  of  ( 6  )  *  that  correspond  to  n  »  1/2,  1.  3/2  are  as  follows: 

zz  y  =o 

Case  9.1(b).  Ineteed  of  using  (9.28),  it  ie  simpler  in  the  caee  of  n«& 
to  compute  (0  )*  directly  from  (9.1) 

ZZ  f)  =0  7  I 

f *'[<*-*>  f, 

- 

These  integrals  are  elementery  giving 


6-JuT*  J<r 


(9.29) 
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which  checks  with  thet  obatined  by  Sneddon  [9j  • 

Case  9.2(b) .  Letting  n  =  1  in  (9.28)  end  carrying  out  the  elgebra 


render 


31) 


which  ie  0(l/r)  for  lerge  r.  At  the  boundery  points  of  r  «=  e,  (€>zz)^  mg  ie 
unbounded. 

Case  9.3(b) .  The  limiting  farm  of  (9.28)  for  n  =  3/2  corresponds 
precisely  to  (8.6)  which  hes  already  been  diecussed. 


10.  Conclusion.  The  problem  of  finding  streeeee  and  displecements  in 
en  elestic  eolid  with  geometric  diecontinuities  hee  been  reduced  to  the 
classical  boundary  problem  of  potentiel  theory.  Hermonic  functions  are 

r 

developed  for  the  case  of  en  elliptically-eheped  plane  of  discontinuity 
whose  faces  are  subjected  to  pressures  end/or  dieplecemente.  While  the  present 
peper  ie  primarily  concerned  with  the  celculetion  of  the  normal  displecement 
w  and  the  corresponding  pressure  ^zz  ,  the  remeining  dieplacements  and 
streeses  mey  be  obteined  from  (2.1),  (2.2),  (2.4),  end  (2.5)  without  difficulty. 
Furthermore,  with  the  knowledge  of  f(x,y,z),  the  stress  components  [2] 
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(10.1) 


Z a  [(t-ZJ)  +  i  ]  , 


(10.2) 

(10.3) 


are  also  known.  The  method  of  solution  outlined  in  the  paper  may  be  used  to 
solve  other  boundary  problems  of  fundamental  interest.  For  example,  the 
skew-symmetric  problem  of  specifying  displacements  u,  v  on  the  surfaces  of 
the  geometric  discontinuity  may  be  formulated  and  solved  in  a  similar  manner. 
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